clear all
version 15
cd "___________"
set more off

scalar idval=7   // Beijing
numlist "0(5)25"
local foo "`r(numlist)'"
do PM10CMLE_withinference
graph combine distest propFx, graphregion(color(white)) xsize(3) ysize(1.5)  rows(1)
graph export "./outputgraphs/main/Figure4B.pdf", replace
do PM10CMLE_robustness
graph combine propFx robust, graphregion(color(white)) xsize(3) ysize(1.5)  rows(1)
graph export "./outputgraphs/main/Figurerobust_Beijing.pdf", replace

scalar idval=67   // Shanghai
numlist "-5(5)5"
local foo "`r(numlist)'"
do PM10CMLE_withinference
graph combine distest propFx, graphregion(color(white)) xsize(3) ysize(1.5)  rows(1)
graph export "./outputgraphs/main/Figure7_Shanghai.pdf", replace
do PM10CMLE_robustness
graph combine propFx robust, graphregion(color(white)) xsize(3) ysize(1.5)  rows(1)
graph export "./outputgraphs/main/Figurerobust_Shanghai.pdf", replace

scalar idval=79   // Tianjin
numlist "-5(5)15"
local foo "`r(numlist)'"
do PM10CMLE_withinference
graph combine distest propFx, graphregion(color(white)) xsize(3) ysize(1.5)  rows(1)
graph export "./outputgraphs/main/Figure7_Tianjin.pdf", replace
do PM10CMLE_robustness
graph combine propFx robust, graphregion(color(white)) xsize(3) ysize(1.5)  rows(1)
graph export "./outputgraphs/main/Figurerobust_Tianjin.pdf", replace

scalar idval=16   // Chongqing
numlist "-5(5)20"
local foo "`r(numlist)'"
do PM10CMLE_withinference
graph combine distest propFx, graphregion(color(white)) xsize(3) ysize(1.5)  rows(1)
graph export "./outputgraphs/main/Figure7_Chongqing.pdf", replace
do PM10CMLE_robustness
graph combine propFx robust, graphregion(color(white)) xsize(3) ysize(1.5)  rows(1)
graph export "./outputgraphs/main/Figurerobust_Chongqing.pdf", replace

